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ABSTRACT 

We calculate the galaxy bispectrum in both real and redshift space adopting the most 
common prescriptions for local Eulerian biasing and Lagrangian evolving-bias model. 
We show that the two biasing schemes make measurably different predictions for these 
clustering statistics. The Eulerian prescription implies that the galaxy distribution 
depends only on the present-day local mass distribution, while its Lagrangian coun- 
terpart relates the current galaxy distribution to the mass distribution at an earlier 
epoch when galaxies first formed. Detailed measurement of the galaxy bispectrum (of 
its reduced amplitude) can help establish whether galaxy positions are determined by 
the current mass distribution or an earlier mass distribution. 

Key words: galaxies: statistics - large-scale structure of Universe. 



1 INTRODUCTION 

Galaxy clustering in the nearby Universe has been mapped 
through a variety of surveys, including different populations 
of luminous objects. These run from optical galaxies in the 
APM, CfA and LCRS surveys to the sources of the IRAS 
catalogue at 60 fim. Reconstructing the overall mass power 
spectrum from these data represents one of the main goals of 
modern cosmology. It is already known, however, that differ- 
ent tracer populations show different clustering amplitudes 
even after redshift-space and small-scale corrections are ap- 
plied. Thus, their clustering patterns are not unambiguosly 
related to any one given mass power spectrum (see Peacock 
1999 for a review). 

The simplest and most common description of biasing 
adopted in the literature is that, at any spatial position x, 
the fluctuation in the number density of galaxies <5g(x) re- 
sponds linearly and locally to the underlying mass fluctu- 
ation <5(x), namely <5 s (x) = & B <5(x), where b E is a space- 
independent bias factor (e.g. Dekel & Rees 1987). As dis- 
cussed below, higher-order bias factors can be introduced, 
but the point is that such a bias prescription is inherently 
Eulerian: it relates the present-day galaxy and mass clus- 
tering properties, ignoring their past evolution. However, if 
gravity is the main force acting in the Universe, there is no 
doubt that galaxy biasing evolves in time, as collapsing mass 
fluctuations keep accreting luminous matter onto them, the 
galaxy distribution eventually relaxing to the mass one (Fry 
1996; Tegmark and Peebles 1998). So, the biasing in the 
present-day galaxy distribution might well be rooted into the 
deep past of the history of the Universe: the strong Lyman 
break galaxy clustering seems to suggest that this might be 
the case (Steidel et al. 1997). Any primordial biasing, aris- 



ing at the epoch of galaxy formation, cannot be described 
by an Eulerian model. Instead, a Lagrangian one has to be 
adopted: it is the primordial fluctuation in galaxies that is 
proportional to the mass fluctuation, S g (q) = b L <5(q), where 
q denotes the Lagrangian position; in general b E differs from 
b L and, in principle, higher order factors can be defined. 

In this Letter we show that the local Eulerian and La- 
grangian bias models are inconsistent. In fact, the cluster- 
ing patterns predicted by the two bias models are different. 
Specifically we study the galaxy bispectrum and skewness, 
both in real and redshift space, on scales where the mildly 
non-linear approximation suffices, starting from Gaussian 
initial conditions. In Section 2 we review the general Eule- 
rian and Lagrangian bias models in terms of infinite hierar- 
chies of bias factors {bf} and {bf}. In the whole discussion, 
these must be considered as free, position-independent, pa- 
rameters. In Section 3 we discuss the galaxy bispectrum and 
skewness in real space, for both bias models. In Section 4 
we carry out the same analysis, but taking into account the 
effect of redshift distortions. Section 5 contains our conclu- 



2 THE TWOFOLD BIASING PRESCRIPTION 

Let us start by fixing the notation of basic quantities. If 
<^ (q) is the primordial gravitational potential (growing 
mode only, smoothed on some scale R and linearly extrap- 
olated to the present time), then 5^(q,z) = D(z) V^o(q) 
is the linear density field, and D(z) its growth factor with 
z the cosmological redshift [we put D(0) — 1]. The linear 
peculiar velocity is given by u^(q) = — V 9 ipo(q), and it 
is constant in time. The Eulerian density field will be indi- 
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cated by <5(x, z) and the n-th order perturbative solutions 
S (n) are such that S = J2„ s(n) (Goroff et al. 1986). The 
Fourier transform is, e.g., S(k) = Jdx.5(x.) expik • x. 



2.1 Local Eulerian bias 

In this approach, the galaxy number density field at a given 
position x and time z (e.g. 'here' and 'now') is assumed to 
be a local function of the underlying mass density field at 
the same location and instant, <5 9 (x, z;R) = £[S(x, z; R)], 
where the smoothing scale 7? is much larger than the typical 
size of the selected objects. Usually, assuming that £ [8] can 
be expanded about 8 — as a power series, an infinite set of 



"Eulerian bias factors" 
1993): 



3=0 J ' 



bj can be defined (Fry & Gaztafiaga 



(1) 



This series is such that (S g ) — and S g (S = —1) = — 1. The 
linear coefficient bf corresponds to the usual bias factor. 
The origin of this local Eulerian prescription is essentially 
phenomenological, and it is a priori devoid of any insight 
about the dynamics of the clustering. Galaxy clustering is 
analyzed for instance in terms of iV-point correlation func- 
tions _ 2 <5 fl (x n , z)}, and the bias factors are tuned to fit 
the observational data. This is the approach that has been 
implicitly adopted in most of the published literature on 
biasing, at least in its leading approximation. 



2.2 Local Lagrangian bias 

According to this alternative prescription, the sites of galaxy 
formation are identified with specific regions of the primor- 
dial density field. It is then appropriate to define a "primor- 
dial" galaxy density field, 5 a (q), measuring the (smoothed) 
overdensity of galaxies in fieri at the Lagrangian position q 
at a given time z (formally z = oo, i.e. 'there' and 'then') 
which is biased with respect to the primordial (linear) den- 
sity field at the same location and instant, namely 5 9 (q) = 
£Mq)l = EJlo^/j'OeoCq)"; here 6 D (q) = VVo(q) is 
the linear density field extrapolated to the present time. We 
can equivalently write, for similarity with the Eulerian case, 



3=0 J ' 



(2) 



with S^'(q,z) = D(z)e (q) and the Lagrangian factors bf 
are defined accordingly in terms of the original 6„j ■ 

Both galaxies and dark matter flow through Eulerian 
space towards mass concentrations. Therefore, even assum- 
ing that the spots of galaxy formation can be identified in 
Lagrangian space, one has to consider large-scale motions 
in order to compute the statistics of present-day structures. 
One therefore needs to assign a dynamical prescription, be- 
cause evolution changes the original galaxy distribution: this 
is the main difference with respect to the Eulerian bias 
scheme, where no dynamics is taken into account. 

One way is to use galaxies as test particles of the under- 
lying gravitational field (Fry 1996), then the evolved galaxy 



density field at the Eulerian position x and instant z is re- 
lated to the primordial galaxy field and evolved density field 
by the relation (Catelan et al. 1998) 



l + S B (x,z) = [l + 5 9 (q)] [l + 5(x,z)] 



(3) 



We stress the fact that eq.(|3|) is inherently non-local. 
Smoothed regions in Lagrangian space can be mapped to 
Eulerian space through the transformation x = q + S(q, z), 
where S is the displacement vector. In the Zel'dovich (1970) 
approximation, the simplest transformation, S(q, z) = 
D(z)u^ (q). Thus, the resulting 5 g (x,z) and 5(x.,z) are 
not deterministically related. In fact, for any given 5 the 
galaxy field 5 g can assume different values (see Dekel & 
Lahav 1999). This stochastic behaviour is inherent to the 
gravitational instability dynamics. 

The question now is the following: are the predictions 
about the clustering (in terms of standard statistics as the 
correlation functions, for example) as deduced from the local 
Eulerian and Lagrangian bias equivalent? In other terms, do 
there exist two sets of non-trivial Eulerian and Lagrangian 
bias factors, {bf} and {bj}, such that the predictions for 
galaxy clustering are identical? In order to answer to these 
questions, let us analyze the galaxy bispectrum from Gaus- 
sian initial conditions as induced by mildly non-linear den- 
sity evolution. 



3 GALAXY BISPECTRUM 
3.1 Eulerian bias case 

The lowest order contribution to the 
galaxy bispectrum (2tt) 3 Sd (ki + k2 +k 3 ) B g (ki, k.2, k^; z) = 
(<5 g (ki, z)5 g {k.2, z)5 g (ks, z)) comes from the appearance of 
non- negligible second-order fluctuations Sg . From eq. (Q), 
this is 5 g 2) (x,z) = bfS (2) (x,z) + ifcf(5 (1)2 (x,z). (Note that 
this expression does not have zero mean, so an offset term 
should be introduced; however, since we are interested in 
the spectral properties of the galaxy clustering, we will 
ignore it since it contributes only to k = 0.) Defining 
ZA2 = ki • k2/fcifc2, and the second-order growth factor 
E ~ — |f2 -2 / 63 D 2 either in an open Universe with no cosmo- 
logical constant (Bouchet et al. 1992) or E fa -|n _1/140 D 2 
in a Universe with a cosmological constant or quintessence 
(Kamionkowski & Buchalter 1999), we introduce the sym- 
metric kernel 

E 



j: 



i- • 



D 2 







, k 2 \ 






K 








-—) 

D 2 J 



^12,(4) 



and the second-order convolution integral operator 

* 8 m (Fry 1984). Finally, we can simply write 5^ 2 ' 
1^J[ 2) and 



u 9 



(bf4 2) 



+ 2&2 



(5) 



Thus, the galaxy bispectrum is (Matarrese, Verde & Heavens 
1997) 



B, 



2D 4 b? 2 



(bfj 



+ c. t. 



(6) 



where P(k, z) = D(z) 2 P(k) is the mass linear power spec- 
trum. The skewness of the galaxy density field smoothed on 
scale R is therefore (Fry 1994), 
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&f(4- 



' D 2 



+ 3bf 



6f 7 (i?) 



bf 2 

,2 ;ji„ r> „„j J2 



(7) 



where 7 = —din a R / din R and cr^ is the rms density on 
scale R. For a scale-free mass power spectrum P(k) oc 
k n and a top-hat smoothing function, one obtains 7 = 
n + 3 (Bernardeau 1994). We remind the reader that in 
the Einstein-de-Sitter Universe, 4 - 2E/D 2 = 34/7. Mass 
bispectrum can be recovered by setting in these formulae 
bf = 1 and 6f = 0, B% = 2D 4 [J { S 2) P{k 1 )P(k 2 )+c.t.] = B m 
(Fry 1984). The growth of B m is self-similar, i.e. mass parti- 
cles do not move from their initial positions, and the wavec- 
tors k actually correspond to those positions. 



3.2 Lagrangian bias case 

Let us now repeat the previous calculations assuming the 
Lagrangian biasing scheme in eq.(Q). In this case, expand- 
ing eq.(^) up to second-order, we obtain, after Zel'dovich 
transforming the Lagrangian coordinate q to the Eulerian 
one x at z, S g = bg + (1 + bg + bf)S w + 8f \ where, 

5 f) = (i + 6 £ )(5 « - D bf u (1) ■ V5 (1) + {bf + \bf)S {1)2 . (8) 

This expression generalizes the analogous one in Catelan et 
al. (1998) for Press-Schechter dark matter halos, for which 
bf = 0. The Zel'dovich approximation, adopted here, suffices 
to transform from q to x: this explains the presence of the 
inertia term, proportional to the velocity. (We assume the 
scale R large enough so that shell-crossing is absent on scale 
R.) The Fourier transform of eq.(0) is 



5(2) = 
where 



((l + tf)4 a) +Bg° 



(2) 



, L . 1 , L . 1 , L f ki k 2 , 



(9) 



(10) 



describes the effects of Lagrangian biasing during the mildly 
nonlinear regime. So, the galaxy bispectrum is 



Bi 



2D 4 (l + bf + bf) 2 



P{ki) P(fc 2 ) + c. t. 



(i+6^)4 2) +4 2) 



(11) 



and the galaxy skewness turns out to be 

QLfm (l + bf)(4-2-^)+6bf+Sbf-(l + bf+bf) 7 (R) 

Ss {R)= o^bJTbJY 2 - (12) 

We emphasize the fact that it is the term r\bf (fei/fe + 
k 2 jk\)v\ 2 in Bg 2 ', Fourier transform of the inertia term 
— fefu' 1 ' ■ in eq.(|^), which carries the signature of the 

gravitational dynamics, inherently absent in the Eulerian 
description. Clearly, such a signature would reflect into a 
distinctive shape dependence of the galaxy bispectrum, best 
quantified by the 'effective' amplitude Q (see below). 



3.3 Disentangling the two biasing schemes 

3.3.1 Galaxy Bispectrum and Skewness 

The different clustering predictions of the two biasing mod- 
els may be emphasized simply by calculating the differ- 
ence AB g of the bispectra in eq.(p|) and eq.(|ll|) or of the 



skewnesses AS g in eq.(Q) and eq.(^). It can be easily veri- 
fied that no two sets of nontrivial independent bias factors 
{bf} and {bf} can be found such that AB g = = AS g . 
We can explicitly write the final expressions assuming that 
at least the lowest-order bias factors are related, namely 
bf — 1 + &o + b\ . This last relation may be easily derived 
in linear regime, but it shows to be preserved even during 
the mildly non-linear regime (Mo & White 1996; Mo, Jing 



White 1997 for the case b% = 0). Thus AB g = B g — B 



AB g = D 4 (l + b% +b^ 2 { 



< 



(1 



2 X 

V12, 



P(k 1 )P(k 2 )+C. t.j 



= AB 9 (1,2) + AB S (1,3) + AB 9 (2,3) 
Correspondingly, the skewness difference is 



AS 9 = 



3(^-^-2^(1 + ^) 



(13) 



(14) 



{1 + b^+bfY 

Intriguingly, the dependence on the filtering scale cancels 
out, and we are left with a residual difference in the skewness 
of the two bias models that is scale-independent. Expres- 
sions for the Einstein-de-Sitter Universe can be recovered 
by setting 1 + E/D 2 = 4/7. 

The two biasing schemes cannot be distinguished by 
measuring the skewness alone, but they can be distinguished 
from the shape dependence of the bispectrum. The impor- 
tant point is that B s , the new term that arises in the 
Lagrangian-biasing bispectrum [cf. eq. (11)] is linearly in- 
dependent of the mass kernel J^ 2 ' and the constant b 2 /2, 
the two terms that make up the Eulerian-biasing bispec- 
trum [cf. eq. (6)]. Thus, no combination of parameters bf 
and bf can allow the Eulerian-biasing bispectrum to mimic 
the Lagrangian-biasing bispectrum. 

Since the bias of high-redshift populations is > 3 while 
that of local populations are closer to unity, we heuristi- 
cally expect bias evolution of order unity, and thus the 
Lagrangian-bias parameters &q and bf to be quantities of 
order unity (further modeling is required to give a more 
precise answer). It is likely that surveys such as the SDSS 
and 2dF will be able to measure the bispectrum with enough 
precision to distinguish the predictions of Eulerian and La- 
grangian biasing if the bias parameters are of order unity. 
Matarrese, Verde & Heavens (1997) have determined that 
SDSS/2dF data should be able to determine bf and bf to 

roughly a few percent within the context of Eulerian-bias 

(2) 

models; that is, the coefficients of Jg and the constant in 
eq. (6) can be determined to a few percent. To distinguish 

between Eulerian and Lagrangian biasing requires that the 

(2) 

data be fit to an additional term, B s , as well. Although 

we have not revisited the calculation in detail, it seems rea- 

(2) 

sonable that if the coefficients of Ji, and the constant can 

(2) 

be fit to a few percent, then the coefficients of B s can be 
fit with a precision not much poorer. In this case, the data 
can discriminate between Eulerian- and Lagrangian-biasing 
schemes. 



3.3.2 Q- Amplitudes 

Eulerian and Lagrangian biasings predict for the galaxy bis- 
pectrum two different shape dependences, which cannot be 
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obtained from one another by simply tuning the bias param- 
eters bf , bf and bf , bf ■ An efficacious way of emphasizing 
such a shape dependence is through the bispectrum 'ampli- 
tudes' Qg ^g , 



and Qg , where, for example, 



Q 3 



£^(ki,k 2 ,k 3 ,z) 



[Pg L (k 1 ,z)Pg L (k 2 ,z) +C.t.] 



(15) 



and similarly for Q g (Fry 1984). Q-amplitudes are essen- 
tially insensitive to the scale and the overall geometry. We 
obtain 

Qf = Q s L + ^r. (16) 

It is very useful to express Qf and Qg in terms of the bispec- 
trum amplitude of the underlying mass density distribution 
Q m = Bm/ ^2 W e have, respectively, 



If 



+ 



bf- 



(17) 



for the Eulerian amplitude (Fry 1994), and, for the La- 
grangian amplitude, 



Qg = 



+ 



Qrt 



+ 



! + + ' (1 
(l + E/D 2 )bf [(1 

( 1 + b L +b L )2 



bf^ + bff 

-^l 2 2)f(fcl)P(fc 2 )+C.t. 

[P(ki)P(k 2 ) + at.] 



(18) 



The novelty contained in eq.(y_8|) is the angular dependence 
appearing in the right hand side: it is not related to the Q m - 
shape dependence, and it is independent from the values 
of {bf}. Measurements of Q g from galaxy catalogs for two 
different shapes of the triangle ki + k 2 + k 3 = can in 
principle disentangle the two biasing factors b\ and 6 2 and 
the two biasing schemes as well. In Figure 1 we plot Qf and 
Qg for a ACDM model; the values of the biasing factors 
bf t 2 and the choice of the scales are based on Scoccimarro 
et'al. (2000). 

More sophisticated and predictive relations may be pro- 
posed if one assumes that the set of Lagrangian bias fac- 
tors {bf} are not free parameters, as in the present discus- 
sion, but rigourously computed within the framework of a 
given theoretical model. The 'excursion set' formalism (Pea- 
cock & Heavens 1990; Bond et al. 1991), for example, where 
dark-matter halos are identified by first-upcrossings of a col- 
lapse thereshold, predicts that bf and bf are functions of 



both halo size and redshift (Mo 
al. 1998). 



White 1996; Porciani et 



4 REDSHIFT DISTORTION EFFECTS 

Given that the two biasing schemes are in principle dis- 
tinguishable, we proceed to calculate the Eulerian and La- 
grangian bispectra in redshift space, which is where they 
are most likely to be measured. Peculiar motions associated 
with structures on any scale distort the clustering pattern in 
redshift space (Kaiser 1987). So, in order to reconstruct the 
actual distribution of galaxies from redshift catalogues, we 
must be able to invert the distortion process. This can be 
easily done if we consider a distant region of the Universe 
so that the distortions essentially occur along the line-of- 
sight, and we restrict to large scales for which the mildly 




Figure 1. The halo bispectrum amplitude Q for configurations 
with sides k\ = 0.05/i6s/Mpc and fc 2 = O.lfoes/Mpc separated 
by an angle 8 for a linear A— CDM power spectrum (fi m = 0.3, 
Ca = 0-7, n = 1, erg = 0.9). The prediction of the local Eulerian 
bias model with the bias parameters estimated from the IRAS 
QDOT 2 Jy redshift catalogue (bf = 0.76, bf = -0.33; Scocci- 
marro et al. 2000) is represented by a dashed line. The continuos 
and dotted lines represents two local Lagrangian bias models with 



: and bf 



1 



-0.24. The value of bf has been fixed 



to match the Eulerian bias prediction for Q at its minimum and 
maximum value (respectively, bf = —0.39 and bf = —0.58): for 
these specific configuration and choice of parameters, discrepan- 
cies between the predictions of the two biasing schemes are about 
10% for the Q-tails and about 50% for the Q-trough. 



non-linear approximation suffices. If r is the physical coor- 
dinate, and u = v ■ r/r is the line-of-sight component of the 
peculiar velocity v, assuming that the observer's peculiar 
velocity is zero, the apparent galaxy fluctuation Sg(s) at the 
apparent position s = (l + u/r)r is related to the actual one 
5 g computed at the same apparent position by the relation 



^(s) = S ai s ) - u '( s ) ~ M ( s )( 5 s( s ) ~ u '( s )) 



(19) 



Here u' indicates the first radial derivative of u. Since in 
this section we will compute the effects of redshift distor- 
tions on the galaxy bispectrum, both for a local Eulerian 
and Lagrangian bias, only corrections up-to second order 
are considered. We remind the reader that in the distant- 
observer limit, the Fourier transform of d/dr —> ik/i where 
fi — k-r/fcr and u = ifj,f(Q)fj/k, where /(fi) ~ O. ' 6 , a is the 
universal scale factor, and rj is the divergence of the velocity 
field. 



4.1 EB galaxy bispectrum in redshift space 

In this case, inserting eq.(^) into eq.(^|), the Fourier trans- 
form of 5g(s) is given by 



(6f +// })5 {1) +1^S 



r(2) cC 2 ) 



(20) 
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where the redshift-distorted symmetric Eulerian-bias kernel 
is 



+ 2 b2 



+ 



2 bl f 



Mi + M2 + Mi A*2 



(zL _i_ zL 

\k 2 ki 



+ f 



2 2 1 / 2 kl 2^2 

Ml M2 + M2^Ml^ + M2^- 



(21) 



(2) 

The quantity ifg describes the second-order contribution 
to rj (Goroff et al. 1986). The distorted galaxy bispectrum 
is (Heavens, Matarrese & Verde 1998) 

B 3 g {E) = 2D 4 (bf + /*?/) (bf + v 2 2f)S ( £ ) P(k 1 )P(k2)+c.t.(22) 



4.2 LB galaxy bispectrum in redshift space 

obtaining, 



We adopt in this case the expression in eq.i 
after analogous calculations, 

S S B = (1 + bo + bf + m 2 /)5 (1) + 1 {2) S {2) . 
Thus, the galaxy bispectrum is 



(23) 



B' g (L) = 2 J D 4 (l + 6,t + 6i +tHf){l + bS +bt + Mi/) 

x 5[ 2) P(fci)P(fc 2 ) +c. t. , (24) 

where, in this bias prescription, the redshift-distorted 
second-order Lagrangian-bias kernel <S^ is 



+ ^u + ^ + bf)/ 



■'2 / ™1 

Ml + M2 + Ml M2 



/fel fc2\ 

U 2 ki) 



+ f 



2 2 

Mi M2 



1 / 2 fel 2 

+ 2MiM2(Mi^+M2^ 



(25) 



4.3 Comparing bias in redshift space 

If we assume, once again, the validity of the algebric rela- 
tion bf = 1 + 6n + bf, it follows that between the redshift- 

(2) C2) T (2} 

distorted kernels holds the relation = + &i J$ + 
ifof — S^?'. This relation should be immediately compared 
with the one in eq.(|l^), to understand that we obtain the 
concise expression between the quantities ABg and AB g 
which emphasize the inconsistency between the two biasing 
prescriptions, 

AB° = (1 + i£ff) (1 + AB g (l, 2) + c. t. , (26) 

where /3 = f /(1 + bo +bj). Thus, the only redshift effect on 
the quantity AB g (i,j) comes from the first-order distortion 
of the galaxy number density field, — (1 + /i 2 (3)8^. It 
has to be like that, if one thinks that the distortion effects 
due to peculiar motions are either independent of the bias 
factors or proportional to the first-order bias factors, then 
they cancel out. In redshift space, eq. (|lrj) becomes 

A R s 

Ql(E)=Ql(L) + ^-^. (27) 

Though the structure of the expression ( |l6| ) is preserved in 
redshift space, Q 9 (Q m )-relations like those in eqs.(|l7|) and 
(^) are not. A comprehensive investigation of the effects of 
redshift distortions on Bf and B g is in progress; see also 
Scoccimarro et al. (1999) for an analysis of B g (E). 



5 DISCUSSION AND CONCLUSIONS 

We compared the galaxy clustering predictions of the local 
Eulerian bias scheme versus those of the Lagrangian one. 
We showed that the two bias models are inconsistent, since 
the predicted three-point galaxy correlations are different. 
A similar inconsistency certainly characterizes correlations 
of higher order, or of lower order but higher perturbative 
corrections. Qualitatively, these results are independent on 
whether the Lagrangian zero-order bias factor b$ is zero, 
as for Press-Schechter dark matter halos, or not, as in the 
most general case we have considered here. The galaxy bis- 
pectrum is much better suited to distinguish between the 
two bias models than the corresponding skewness, since the 
latter is spatially averaged: the bispectrum depends on the 
shape of the triangle ki + k2 + k3 = 0, thus two shapes 
can disentagle the two bias factors b± and 62 (Matarrese, 
Verde & Heavens 1997; Scoccimarro 2000) and the two bias 
models. The shape dependence is best quantified by the Q- 
amplitudes discussed in Subsection 3.3.: the reduced am- 
plitude of the bispectrum from Lagrangian bias in eq.([l8|) 
displays a dependence on the triangle configuration which 
is not contained in eq.(^), and, since it is dynamically in- 
duced, which is independent from the bias factors. The next 
generation redshift catalogues, as the ongoing Two Degree 
Field Survey and the Sloan Digital Sky Survey, will contain 
enough galaxies to establish whether B g {E) in eq.(p2j) or 
B g (L) in eq.(p4|) better fits the observational data, but they 
cannot be both correct, whatever the assumed cosmology. 

Both bias schemes represent rather extreme and ideal- 
ized approaches. Lagrangian models imply a sort of infinite- 
memory process, since the sites for galaxy formation are 
known from the beginning, and dynamical evolution changes 
their spatial distribution. On the other hand, in local Eule- 
rian schemes galaxies are simply 'painted' on a snapshot 
of the density field, without a record of the past. How- 
ever, even though real galaxy formation is probably a pro- 
cess with intermediate characteristics with respect to the 
biasing schemes discussed here, recent models based on a 
Lagrangian selection of the sites for object formation were 
shown to be very successful in reproducing the clustering of 
dark-matter halos found in numerical simulations (e.g. Cate- 
lan, Matarrese & Porciani 1998; Porciani, Catelan & Lacey 
1999). The issue discussed in this Letter surely deserves fur- 
ther investigation, both in real and redshift space. It would 
be of interest to test which biasing scheme better describes 
galaxy power spectrum and bispectrum from a combination 
of numerical simulations and semianalytic models (Porciani 
et al., in preparation). 
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